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We propose a scheme to simulate and explore Weyl semimetal physics with ultracold fermionic 
atoms in a two-dimensional square optical lattice subjected to experimentally realizable spin-orbit 
coupling and an artificial dimension from an external parameter space, which may increase experi¬ 
mental feasibility compared with the cases in three-dimensional optical lattices. It is shown that this 
system with a tight-binding model is able to describe essentially three-dimensional Weyl semimetals 
with tunable Weyl points. The relevant topological properties are also addressed by means of the 
Chern number and the gapless edge states. Furthermore, we illustrate that the mimicked Weyl 
points can be experimentally detected by measuring the atomic transfer fractions in a Bloch-Zener 
oscillation, and the characteristic topological invariant can be measured with the particle pumping 
approach. 
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I. INTRODUCTION 

Weyl semimetal (WSM), as an exotic topologically 
nontrivial state of matter in three dimensions (3D) [IH8], 
has a singly-degenerate band structure with paired bulk 
band crossings, named as Weyl points. A linear disper¬ 
sion relation in all three momentum directions exhibits 
near a Weyl point, with low-energy excitations resem¬ 
bling the well-known Weyl fermions in particle physics. 
Each Weyl point has a topological charge, namely the 
Chern number on the gapped 2D sphere enclosing a Weyl 
point in momentum space, and a pair of points with 
opposite charges support topologically protected gapless 
Fermi arc surface states ®m- The Weyl points and the 
Fermi arcs in a WSM are fundamentally interesting, and 
are expected to give rise to exotic phenomena absent in 
fully gapped topological phases, such as anomalous elec¬ 
tromagnetic responses 0 ED Eg. Recently, more and 
more efforts have been made on experimental realization 
of WSMs in real materials HM3 and in artificial sys¬ 
tems, such as analogous Weyl points in photonic crystals 

mm- 

On the other hand, ultracold atomic gases in optical 
lattices provide a powerful platform to simulate various 
quantum states of matter (HI EE], In particular, recent 
experimental advances in engineering spin-orbit coupling 
and artificial gauge field for ultracold atoms [22H261 have 
pushed this system to the forefront for studying topolog¬ 
ical states of matter HMD- Another recently developed 
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technique related to the realization of topological phases 
in cold atom systems (also in photonic quasicrystals, e.g. 
see Ref. m) is extending an artificial dimension in op¬ 
tical lattices provided by an external cyclical parameter 
[43l 44] or the internal atomic degrees of freedom [05 , 46 . 
Using this technique, one can study the physics of topo¬ 
logical states in optical lattices attributed to dimensions 
higher than their own. Nowadays, the Zak phase in topo¬ 
logically nontrivial Bloch bands realized in ID optical lat¬ 
tices has been measured [39]. The 2D Chern insulators 
and the Hofstadter bands have been also realized exper¬ 
imentally in 2D optical lattices with the Chern number 
therein being successfully probed @DJ E]. The exper¬ 
imental observation of chiral edge states in ID optical 
lattices subjected a synthetic magnetic field and an ar¬ 
tificial dimension has been reported [46] . By stacking 
multilayers of 2D atomic Chern or Hofstadter insulators, 
it was proposed to realize WSMs in 3D optical lattices 
[47U49] . This construction method was extended to simu¬ 
late WSMs in ID double-well optical lattices with two de¬ 
grees of freedoms modulating the hopping terms [50] , but 
it seems extremely hard to independently control them 
in practical experiments. Due to the elusive nature of 
WSMs and their intrinsic exotic properties, other feasible 
schemes for their experimental realization or simulation 
would be still of great value. In particular, the simula¬ 
tion and detection of the Weyl points and the relevant 
topological properties of WSMs with an artificial dimen¬ 
sion in feasible 2D cold atomic systems are still badly 
awaited. 

In this paper, we propose a scheme to simulate and ex¬ 
plore essential physics of WSMs with ultracold fermionic 
atoms in a 2D square optical lattice, subjected to ex¬ 
perimentally realizable spin-orbit coupling and an artifi¬ 
cial dimension from an external cyclical parameter. We 
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first show that this system with a tight-binding model 
is able to mimic 3D WSMs with tunable Weyl points. 
The topological properties in this system are then fur¬ 
ther investigated by calculating the out-of-plane momen¬ 
tum, say k z , -dependent Chern number and the gapless 
edge states. Finally, we propose practical methods for 
the experimental detection of the mimicked Weyl points 
and the characteristic topological invariant in the pro¬ 
posed cold atomic system. With numerical simulations, 
we demonstrate that the analogous Weyl points in this 
system can be sharply probed by measuring the Landau- 
Zener tunneling to the excited band after a Bloch oscil¬ 
lation and the k z -dependent Chern number can directly 
be extracted from the center shift of the hybrid Wannier 
functions, both of which are measurable with cold atoms 
in optical lattices by the time-of-flight images. Compar¬ 
ing with the schemes evolving complex synthetic mag¬ 
netic fluxes in 3D optical lattices [T7H49] . our scenario 
not only is different but also may increase experimental 
feasibility. In addition, the proposed system with high 
tunability and convenient detection ways may provide a 
promising platform for exploring exotic WSM physics. 

The paper is organized as follows. Section II intro¬ 
duces the cold atom system with a tight-binding model 
for simulating the WSM phase with tunable Weyl points. 
In Sec. Ill, we elaborate the topological properties in 
this model by calculating the /^-dependent Chern num¬ 
ber and gapless edge states. In section IV, we propose 
practical methods for the experimental detection of the 
mimicked Weyl points and the characteristic topological 
invariant in the proposed cold atomic system with nu¬ 
merical simulations. Finally, a short conclusion is given 
in Sec. V. 


II. MODEL 

In this section, we construct a tight-binding model in a 
2D square lattice subjected to an artificial dimension for 
simulating the tunable WSM phase. Let us consider a 
non-interacting spin-1/2 ultracold degenerate fermionic 
gas (labeled as spins t and loaded in a 2D square op¬ 
tical lattice in the xz plane. In the tight-binding regime, 
the atomic hopping between two nearest-neighbor lattice 
sites can be spin-conserved hopping or spin-flip hopping 
which can be achieved by Raman coupling between the 
two spin states | t) and | i) [26] . The tight-binding 
Hamiltonian of this cold atom system is considered to be 

n = H L + Hsoc + H P , ( 1 ) 

which consists of the spin-conserved hopping in the 2D 
lattice Hl, the spin-flip hopping along the x axis Hsoc , 
and Hp representing additional external coupling poten¬ 


tial: 

Hl = - X t a (at CT a i+ jj )CT + h.c.) , (2) 

i,cr ,?7 

H-SOC = ts ^ ' (aj[^• < -'’) ?(^) 

H P = T Z J2 («i, t ai -t - «u®u) ( 4 ) 

i 

+ X (“i,t“U + “u^t) • 

i 

Here di ?cr (dj a ) annihilates (creates) a fermion on site 
i = (i x ,i z ) with spin a = {t, 1} and r] = { x , z}. The spin- 
dependent hopping amplitude t a in Hp is assumed to be 
1 1 = —t± = to, which can be achieved in experiments 
by Raman laser or modulation engineering hopping m, 
such as laser-induced Peierls phase factor e l7r for the 
| 1) atomic hopping. The spin-flip hopping term Hsoc 
can be induced by the experimentally realized spin-orbit 
coupling with equal Rashba and Dresselhaus amplitudes 
through a two-photon Raman process [23H26] . Finally, 
the Hamiltonian Hp is composed of a Zeeman potential 
and an additional coupling term on the spin states. The 
Zeeman potential is generated by a external magnetic 
field and the coupling term is induced by additional Ra¬ 
man lasers m- The Zeeman and Raman terms can be 
turned fully and independently in cold atom experiments 
with long coherence time m , and thus we choose the pa¬ 
rameters in the following forms 

T z = m z - 2t 0 cos<9, T x = 2t s sm0 , (5) 

where m z is a constant parameter and 0 is a cyclical 
parameter that can vary from 0 = —i r to 0 = 7r. 

Under the period boundary condition, the total Hamil¬ 
tonian of the system can be rewritten as 

^=X«L mk)l a ,a ka ,, (6) 

k,cr<j' 

where a^ G = 1/VV e ~ tk is the annihilation op¬ 
erators in momentum space k = ( k x ,k z ), and H(k ) is 
the Bloch Hamiltonian given by H(k ) = d(k ) • <r. Here 
a = ((T x ,cry,<j z ) are the Pauli matrices acting on the 
atomic components and d = (d x ,d y ,d z ) is the Bolch 
vectors: d x = 2£ s sin$, d y = 2£ s sin£r c , and d z = 
m z — 2to cos k z — 2to cos k x — 2to cos #, with lattice spacing 
a = 1 and h = 1 hereafter. The energy spectrum of the 
system is then given by E±(k,0) = ±|d(fe,0)|. Treat¬ 
ing the cyclical parameter 0 as the Bloch momentum k y , 
which acts as an artificial dimension of the y axis under 
periodical condition along this direction, we can consider 
the 2D system in a 3D parameter space with the 3D- 
extended momentum k = (k x ,0,k z ). The bulk energy 
bands are fully gapped except the points that satisfy the 
following conditions: 

sin k x = sin 0 = 0, (7) 

m z — 2£q(cos k z + cos k x + cos 0) = 0, (8) 
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which determine the boundaries between the mimicked 
WSM phase and the bulk insulating phases as discussed 
in the following. 

We now show that the WSM phase can be simulated in 
this system. Without loss of generality, we focus on the 
parameters to > 0, t s > 0 and m z > 0 in the system in 
the rest of the this paper. The first condition of Eq. (7) 
requires k x = {0,7r} and 0 = {0,7r}, and consequently 
the second condition of Eq. (8) reduces to: 

m z — 2to cos k z = 0, for (k x , 0) = (0, i r) and ( 7 r, 0); (9) 

m z — 2to cosk z — 4to = 0, for (^,6^) = (0,0); (10) 

m z — 2to cos k z + 4to = 0, for (k x , 0) = ( tt , 7 t ). (11) 

When one of the three equations (9-11) is fulfilled with 

the double-valued k z G [—7r,7r], a pair of degenerate 
points will exhibit in extended 3D band structure, which 
can be regarded as analogous Weyl points in the k space. 
For the considered parameter regime, Eq. (11) can not 
be fulfilled, while Eq. (9) gives rise to two pairs of analo¬ 
gous Weyl points denoted by W = (0, 7 r, db arccos ^j) 
and W 2 ,± = ( 77 ,0, ± arccos |^-) for m z < 2£o- Eq. 
(10) gives rise to a single pair of Weyl points W± = 
(0,0, ± arccos mz 2 ~ Q 4t ° ) for 2to ^ m z ^ 6£o- In Fig. 1, we 
plot the extended 3D bulk band structure E + {k x , #, k z ) 
for the typical parameters m z with fixed hopping energy 
to = t s = 1. Therefore, the mimicked Weyl points can 
be found in the positions where = 0. When m z = 1 
[Fig. 1(a)], two pairs of analogous Weyl points locate at 
(k x ,0,k z ) = (0,7r, dz7r/3) and ( 77 ,0, ±7r/3) as expected, 
and the critical case appears when m z = 2 [Fig. 1(b)], 
with the two pairs merging at (0,7r,0) and ( 77 ,0,0), and 
another pair being created at (0,0, ±7r). When m z = 4 
[Fig. 1(c)], the newly created pair of analogous Weyl 
points move to the position (0, 0, ± 7 t / 2 ), and then merge 
at the location (0, 0, 0) when increasing the parameter 
to the critical value m z = 6. From Fig. 1, we can 
find that the number and location of the analogous Weyl 
points are able to be tuned just by varying the parameter 
m z . This demonstrates that one can conveniently create, 
move, and merge the mimicked Weyl points in this cold 
atom system. 

Without loss of generality, hereafter we focus on the 
simplest case of a singe pair of Weyl points W ±, where 
the bulk bands only touch at the two distinct points 
in the k space. We can expand the Bloch Hamiltonian 
around the two mimicked Weyl points up to liner order in 
k , and then obtain the low-energy effective Hamiltonian 

,± = V x Qx&y T VyQyCf x A VzQz&zi (12) 

where v x = v y = 2t s and v z = 2to are the effective Fermi 
velocity, and q = (q x ,q y ,qz) = k — W± for the two 
points, respectively. The Hamiltonian (12) is an analo¬ 
gous anisotropic Weyl Hamiltonian for the Weyl fermions 
and can be written as Hw,± = j Qi a ij a ji where [a^-] is 
a 3x3 matrix with elements a xy = a yx = 2£ s , a zz = ±2£o 
and zero otherwise. Thus the chirality of the two Weyl 



FIG. 1: (Color online) Color-coded extended 3D bulk band 
structure E+(k x , 6, k z ) for the typical parameters m z . (a) 
m z = 1 with two pairs of analogous Weyl points located at 
{k x ,6,k z ) — (0,7t,±7t/3) and ( 77 , 0, ±7t/3); (b) m z = 2, criti¬ 
cal case with two pairs of analogous Weyl points merging at 
(0, 7T, 0) and (7r, 0, 0), and another pair appear at (0,0, ±tt); 
(c) m z — 4 with one pair of analogous Weyl points located at 
(0,0,±7t/2) ; and (d) m z = 6, critical case with the pair of 
analogous Weyl points merging at (0, 0, 0). Other parameters 
in (a-d) are to = 1 as the energy unit and t s = 1. 


points W± can be defined as x± = sign(det[a^-]) = ±1, 
respectively. 


III. TOPOLOGICAL PROPERTIES 

The opposite chirality of the paired analogous Weyl 
points reflects the topological properties of this system, 
which can be characterized by the corresponding topolog¬ 
ical charge. Actually, a mimicked Weyl point here can be 
regarded as a monopole in the k space, whose topological 
charge equals its chirality mm- The paired Weyl points 
are robust to perturbations which add a cr x , y , z term to 
the Hamiltonian because their charge is quantized: the 
monopole charge cannot vary under a continuous change 
of the Hamiltonian and they can only change at creation 
or annihilation of a monopole-antimonopole pair. More¬ 
over, it was shown recently that the gapless edge modes 
of a WSM are also robust against disorders |8j. 

To further look into the topological properties of this 
system, we consider the Hamiltonian (1) with the so- 
called dimension reduction method. We treat k z as an 
effective parameter and reduce the original system to a 
(/^-modified) collection of effective 2D systems, as tight- 
binding chains along the x axis. Such a reduction method 
is valid for the bulk system and the reduced chains if k z 
is a good quantum number. For a fixed k z , the reduced 
Bloch Hamiltonian Et(k x ,6,k z ) —)► Elk z (k x ,6), with 

H kz ( k x , 9) = 2 t s sin 0a x + 2 t s sin k x a y (13) 

+ [M z — 2to cos k x — 2to cos 0)a z 
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acting as an effective 2D band structure in the k x — 9 
parameter space, where M z = m z — 2to cos k z . Since the 
extended 3D bulk bands are fully gapped when k z ^ ±k z 
with k c z = arccos[(ra z — 4£o)/2£o], in this case, Hk z (k x ,0) 
describes a system of the effective 2D Chern insulator. 
Therefore, one can have a well-defined Chern number 

C kz = T J dk x J do d- (d kx d x dgct) 

_ ( —sign (M z ), 0 < \M Z \ < 4 1 0 ; , , 

“\0, \M Z \>U 0 . ^ i4j 

where d = d/\d\. For our case within the parameter 
regime 2 £q ^ rn z ^ 6to, we can obtain that Ck x = — 1 
when k z G (—k z ,k z ) and otherwise Ck z = 0. For each 
k z -fixed chain with a nonzero Chern number, according 
to the bulk-edge correspondence and the index theorem 
[52], there must be topologically protected gapless if an 
edge is created, i.e. finite lattice sites along the x axis. 
This reduced system of 2D analogous Chern insulator is 
described by the reduced Bloch Hamiltonian Hk z (k x ,6), 
and the corresponding tight-binding chain Hamiltonian 
is given by 

fie = - ( a i.+i,t a< -.t - «t+1,^*4)+ h - c - 

i x 

^ ^ t's —1,4^ A h.c. 

ix 

+ Y2( M z~ 2t 0 cos 0) (at >t ai x , t - 4*4®** ,4-) 

^X 

+ J2 2 1* sin6> (al xt d ix 4 + at^a^f) . (15) 

ix 

In Fig. 2, we numerically calculate the energy spec¬ 
trum of the reduced chain with length L x = 40 under 
open boundary conditions for different parameters. From 
Figs. 2(a-c), we can see the variation of the energy spec¬ 
trum by changing the parameter k z for fixed m z = 4to. 
For k z = 0 in Fig. 2(a), the spectrum contains two sym¬ 
metric bands with an energy gap opened, which accom¬ 
panies some in-gap modes splitting into two branches. 
The two branches of edge modes cross and connect the 
separated bands. It indicates that the system with the 
lower band filled is a topologically nontrivial insulator 
with Chern number Ck z = — 1. As we continuously in¬ 
crease the parameter k z: the band gap closes at the crit¬ 
ical value k z = k z = 0.57T [Fig. 2(b)] and then reopens 
without the in-gap edges modes [Fig. 2(c)], indicating 
the system in the topologically trivial case with Chern 
number Ck z = 0. Such a /^-parameter induced transi¬ 
tion between topologically distinct regimes in numerical 
simulations is consist with the previous analytical calcu¬ 
lations. In Figs. 2(d) and 2(e), we show the spectrum of 
topologically nontrivial states for k z = 0 and k z = 0.57T 
with fixed m z = 3£o, in which cases the critical value 
is k z = 27r/3. Despite the shapes of the spectra be¬ 
ing somewhat different, the existence of continuous edge 
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FIG. 2: (Color online) Energy spectrum of the reduced tight- 
binding chain of analogous Chern insulator with lattice sites 
L x = 40 under open boundary conditions, (a) non-trivial 
case for k z — 0 and m z =4; (b) critical case for k z = 0.57T 
and m z = 4; (c) trivial case for k z — tt and mn z — 4; (d) 
non-trivial case for k z — 0 and m z =3; (e) non-trivial case 
for k z — 0.57T and m z = 3. The insert figure in (a) shows 
the density distribution of two typical edge modes. Other 
parameters in (a-e) are to = 1 as the energy unit and t s — 1. 


modes connecting the lower and upper bands indicates 
that the system is topologically nontrivial in these cases. 
The density distribution (profile) of the edge modes in 
these topologically nontrivial cases are similar to the two 
typical ones shown in the insert figure in Fig. 2(a), and 
they will gradually spread into the bulk when their ener¬ 
gies are closer to the bulk bands. 

In the context of 3D WSM phase, a pair of Weyl points 
can be viewed as a monopole-antimonopole pair in the 
momentum space. As a consequence of these monopoles, 
there are chiral Fermi surfaces with the energy E = 0 
connecting the Weyl points on the surface not perpen¬ 
dicular to the z axis, named as Fermi arc jam. In this 
2D system with an artificial dimension, in a similar way 
by assuming that 6 and k z are good quantum numbers, 
we can consider edge modes along the x axis in the mim¬ 
icked WSM phase. The system in this case returns to 
the reduced chain described by Hamiltonian (15) and we 
find that analogous Fermi-arc zero modes emerge by nu¬ 
merically calculating the energy spectrum of the edges 
modes, as shown in Fig. 3. In this case, the pair of 
the analogous Weyl points are connected with straight 
Fermi lines consisting of gapless zero-energy edge modes, 
which are analogous to Fermi arcs in the 0 — k z parameter 
space (denote by green dashed lines). The correspond- 
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FIG. 3: (Color online) Energy spectrum of the edge states in 
the mimicked WSM phase for (a) m z — 4 and (b) m z — 3. 
The green dashed lines denote straight Fermi lines formed 
by gapless zero-energy edge modes, which are analogous to 
Fermi arcs connecting the pair of the mimicked Weyl points 
W±. They are located at ( 6,k z ) = (0, ±7 t/2 ) in (a) and 
(0,k z ) = (0,±27r/3) in (b), respectively. The regimes of left- 
edge and right-edge modes in the parameter space denoted 
by blue ’L’ and red ’R’ are {E < 0, 6 < 0} |J{^ > 0,0 > 0} 
and {E > 0,0 < 0} |J{^ < 0,0 > 0}, respectively. The edge 
modes have similar profiles with those shown in Fig. 2(a). 
Other parameters in (a) and (b) are to = 1 as the energy unit 
and t s = 1. 


ing regimes for left-edge and right-edge modes are also 
shown in Fig. 3(a) and 3(b). The profiles of these edge 
modes are similar with those shown in Fig. 2(a), and the 
edge modes closer to the Weyl points and the bulk bands 
spread more into the bulk than those in the center of the 
parameter regime. 


IV. EXPERIMENTAL DETECTION OF 
MIMICKED WSM PHASE 

So far, we have introduced the system for simulating 
the WSM physics and explored the relevant topological 
properties. In this section, we propose practical meth¬ 
ods for their experimental detection in this cold atom 
system. We first show that the mimicked Weyl points 
can be probed by measuring the atomic Zener tunneling 
to the excited band after a Bloch oscillation, and then 
propose a feasible technique to obtain the k z -dependent 
Chern number directly from the center shift of the hybrid 
Wannier functions. 


A. Detection of the Weyl points 

The band touching points can be monitored from the 
atomic fraction tunnelling to the excited band in Bloch 
oscillations, as recently experimentally demonstrated to 
probe the Dirac points and the topological phase transi¬ 
tion in a honeycomb optical lattice [53, 54. This Bloch- 
Zener-oscillation technique can be extended to detect the 


simulated Weyl points in 3D optical lattices [49] and in 
our system. Let us consider a 2D cloud of noninteract¬ 
ing fermions prepared in the ground band in this sys¬ 
tem and an external constant force F is applied along 
the tj (rj = x,z) direction, which pushes the atoms mov¬ 
ing along the k v direction. For fixed parameters 6 and 
ra 2 , one can observe the quasi-momentum distribution of 
the transfer fraction in the excited band from the time- 
of-flight measurements after a Bloch cycle [551 [54] . One 
can perform such a measurement for varying parameters, 
and finally the mimicked Weyl points in the extended 3D 
Brillouin zone can be verified. 

Without loss of generality, here we consider the case of 
a singe pair of Weyl points described by the low-energy 
effective Hamiltonian (12). In this case, the transfer frac¬ 
tion £ x (k z ) [€z(k x )] along the k x (. k z ) direction for k z - 
trajectory (/^-trajectory) can be written as [551 El] 


Uk Z ) = P X LZ {k Z ), 


Cz (Me) 


2P£z(Me) 


1 — P[z(kx) 5 


(16) 

(17) 


where = e - 7TA l( k z)/^v x F an q pz^ _ e -irA 2 z (k x )/4v z F 
are the Landau-Zener transition probabilities, with = 
2 E + {k x = 0 ,k z ) and A^ = 2 E+(k x ,k z = k c z ) denoting 
the energy gaps for the Landau-Zener events along the 
k v direction. For simplicity, here we treat trajectories for 
different k^ as independent by neglecting the harmonic 
trap and assuming the two tunneling events along the k z 
direction incoherent. Such a simplified model was shown 
to be sufficient for comparison with the experiments un¬ 
der realistic conditions mm- 

In Fig. 4, we calculate the transfer fractions £ XjZ us¬ 
ing Eqs. (16) and (17) for some typical parameters. The 
quasi-momentum distributions £ x (k z ) for varying param¬ 
eter 6 in Figs. 4(a) (with m z = 3) and 4(b) (with 
m z = 4) show two maximum transfer positions in the 
whole k z — 0 space, which correspond to the k z and 6 
positions of the paired Weyl points as excepted. These 
peaks of the maximum transfer factions are sharp since 
the transition probability in a single Landau-Zener event 
increases exponentially as the energy gap decreases. In 
other words, the emergence of the mimicked Weyl points, 
at which the band gap closes, is accompanied by a dra¬ 
matic increase in the transfer fraction £, x (k z ,0). In con¬ 
tract, there are two subsequent Landau-Zener transitions 
along the k z direction. This results in the distribution of 
the final transfer £ z (k x ,Q) exhibiting the ring-type pro¬ 
file of the maximum value £ z = 0.5, as shown in Figs. 
4(c) and 4(d). The position of maximum dip with the 
value £ z ~ 0 inside the ring profile of £ z (k x ,0) indicates 
Me = 0 = 0 for the mimicked Weyl points in these cases. 
We also plot the distribution £ x (k z ,m z ) for fixed 6 = 0 
in Fig. 4(e), where the maximum transfer positions of 
£ x (k z ) correspond well to the expected k z positions of the 
mimicked Weyl points as plotted by the dashed line. The 
ix{k z ,6) distributions for different t s will be modified, 
but the peaks of the transfer fractions remains, with an 
example shown in Fig. 4(f). Therefore, combining with 
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FIG. 4: (Color online) The transfer fractions for tunable pa¬ 
rameters. The quasi-momentum distribution £ x (k z ) for vary¬ 
ing parameter 0 with (a) m z — 3 and (b) m z — 4; the quasi¬ 
momentum distribution £ z (k x ) for varying parameter 6 with 
(c) m z = 3 and (d) m z = 4; (e) £ x (k z ) for different parameter 
m z with 0 = 0; (f) £ x (k z ) for different parameter 6 with m z 
and t s = 0.5. The white dashed line in (e) shows the expected 
k z positions of the mimicked Weyl points W ± for varying m z , 
which are well agree with the maximum transfer positions of 
£,x{kz). The maximum £ x (k z , 0) in (a) and (b) also correspond 
to the k z and 6 positions, as shown in Fig. 3. The position of 
maximum dip inside the ring profile of £ z (k x , 0) in (c) and (d) 
indicates k x = 0 = 0 for the Weyl points. Other parameters 
in (a-f) are to = 1 as the energy unit, F — 1, and t s — 1 
except that t s — 0.5 in (f). 


the measurements of transfer distributions £ x (k z ,6) and 
€z(k x ,0) from Bolch-Zener transitions along the k x and 
k z direction, one can resolve the mimicked Weyl points 
in the k space for different m z . This detection method 
is robust against fluctuations of the hopping strengths, 
which are engineered by laser beams in this cold atom 
system. In addition, since the creation or annihilation of 
the analogous Weyl points by external parameters (such 
as m z ) will accompany topological phase transition, then 
the critical point of the transition can be identified by us¬ 
ing this method. 


B. Detection of the /^-dependent Chern number 

In solid-state materials, the Chern number value can 
be revealed from a routine measurement of the Hall con¬ 
ductance, however, the detection method is usually quite 
different in cold atom systems [30] l33ll36l I391I4T1 [55] . 
Since the Chern number comes from the integral of Berry 
curvature over the first Brillouin zone, one way to re¬ 
veal the topological order is to measure the Berry curva¬ 
ture. It was theoretically proposed and experimentally 
performed in optical lattices to map the Berry curvature 
from the transverse drift induced by the anomalous ve¬ 


locity in Bloch oscillations pi I3H HH1 EH] and from an 
atomic interferometry in momentum space [S3 ESI ESI- 
Besides, it was also proposed to measure the Berry cur¬ 
vature from atomic momentum distributions for different 
spin and spin-mixing components [301 ESI- These meth¬ 
ods for measuring the Berry curvature could be poten¬ 
tially extended to our proposed cold atom system, as one 
can perform similar measurements for a fixed value of the 
parameter 0. 

Below we propose a practical method to directly probe 
the /^-dependent Chern number given by Eq. (14) in our 
system, based on a generalization of topological pumping 
in optical lattices [56H58] . Remind that in order to ob¬ 
tain the topological invariant, we treat our 2D system as 
a collection of tight-binding chains along the x axis de¬ 
scribed by Hamiltonian (15), which can realize the anal¬ 
ogous Chern insulators defined in the k x — 0 space as 
different slices of parameter k z . The polarization of this 
ID insulator can be expressed by means of the centers 
of the hybrid Wannier functions (HWFs) [59], which are 
localized in the x axis retaining Bloch character in the 
6 artificial dimension in our case. This polarization de¬ 
pends on the parameter k z and is a function of #, which 
acts as an external parameter under which the polariza¬ 
tion changes. When 0 is adiabatically changed by 2i r, the 
change in polarization, i.e., the shift of the HWF center, 
is proportional to the Chern number [59l [60] . This is 
a manifestation of topological pumping [ 511 ], with 0 be¬ 
ing the adiabatic pumping parameter. In our system, the 
HWF center for a tight-binding chain in Hamiltonian (15) 
is given by m 


(n x (0,k z )) 


^2i x i'xpij'xi kz) 

5 ~2i x ptyxi kz) 


(18) 


where i x is the lattice-site index, and p(i x ,Q,k z ) is the 
density of the HWF and denotes the atomic densities 
resolved along the x direction as a function of 0 and k z . 
Here the hybrid density can be written as 


pij'x , 0 , k z ) — ^ ^ \^x , @1 k z ) (i x , k z |, ( 19 ) 

occupied states 


where | i x ,0,k z ) is the hybrid eigenstate of the system. 
Experimentally, p(i x ,0,k z ) can be measured by the hy¬ 
brid time-of-flight images Ml referring to an in situ mea¬ 
surement of the density distribution of the atomic cloud 
in the x direction during free expansion along the z di¬ 
rection for a fixed 6. 

We numerically calculate ( n x (6 , k z )) in a tight-binding 
chain of length L x = 100 at half-filling (assuming the 
Fermi energy Ep = 0), with the results for typical pa¬ 
rameter m z being shown in Fig. 5. For the case of m z = 4 
in Figs. 5(a) and 5(b), the expected critical values be¬ 
tween the trivial insulator with Ck z = 0 and Chern insu¬ 
lator with Ck z = — 1 are k z = =b k c z = ±0.57r. As shown 
in Fig. 5(a), (n x (0)) shows discontinuous jumps of one 
unit cell within the region k z G (—0.57T, 0.57t), which dis¬ 
appear outside this region. To be more clearly, we plot 
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FIG. 5: (Color online) The HWF centers in a tight-binding 
chain of length L x — 100 at half-filling as a function of the 
adiabatic pumping parameter 0 for different k z . (a) The pro¬ 
file ( n x (6,k z )) for the parameter m z — 4, where ( n x (0 )) (do 
not) shows a jump of one unit cell for k z (outside) within the 
region (—0.57T, 0.57 t), with typical examples shown in (b). (c) 
The profile for m z = 6.1 shows no jump of (n x (0)) for all k z , 
with typical examples shown in (d). The corresponding /in¬ 
dependent Chern number Ck z is also plotted both in (a) and 
(c), with the white dashed lines denoting the critical value be¬ 
tween the trivial ( Ck z — 0) and nontrivial ( Ck z = —1) cases. 
Other parameters in (a-d) are to = 1 as the energy unit and 

t a = 1 . 



FIG. 6: The maximum variation of the HWF-center shift 
5n x . (a) Sn x as a function of k z for m z = 4. (b) 5n x as a 
function of m z for k z — 0. Other parameters in (a) and (b) 
are L x — 100, to = 1 as the energy unit, and t s — 1. 


{nx(6)) fo r k z = 0 and k z = 0.67T as two examples in Fig. 
5(b). The one-unit-cell jump driven by 6 indicates that 
a single particle is pumped across the system m edi , 
as expected for Ck z = — 1. The results for m z = 6.1 
show no jump of (n x (0)} for all k z in Figs. 5(c) and 
5(d), which is consistent with the expected Ck z = 0 for 
m z > 6 in this case. This establishes a direct and clear 
connection between the shift of the hybrid density cen¬ 
ter and the topological invariant. Thus one can directly 
extract the /^-dependent Chern number from the hybrid 
time-of-flight images in this cold atom system [57]. The 
proposed detection method based on topological pump¬ 
ing is not affected by a weak harmonic trap as it does not 


require the presence of the sharp edge states [571 [58]. 

In the HWF-center shifts (n x (6)) as shown in Fig. 5, 
the one-unit jump for topologically nontrivial cases and 
the continuous variation for trivial cases may be hard to 
distinct in realistic experiments. It will be more practical 
to extract the topological invariant from the maximum 
variation of (n x (6)) during adiabatically varying 0 

5n x = max{(n a ,(6>))} - min{(ra x (0))}. (20) 

In Fig. 6. we numerically calculate Sn x as a function of 
k z ( m z ) with fixed m z = 4 ( k z = 0) for a tight-binding 
chain. One can find that Sn x ~ 1 in the parameter 
regimes where the system is topologically nontrivial, i.e. 
k z E (—0.57T, 0.57r) in Fig. 6(a) and m z < 6 in Fig. 6(b), 
while 5n x ~ 0 outside these regimes, corresponding to the 
trivial cases. Therefore, as long as the system is not close 
to the critical points, such as k z = ±0.57r in Fig. 6(a) and 
m z = 6 in Fig. 6(b), one can extract the topological in¬ 
variant in such a measurement as 5n x = | Ck z | • This sim¬ 
ple relationship does not well survive near the topological 
transition points, where the topological invariant will be 
more difficult to determine in experiments. It is inter¬ 
esting to note that the proposed measurement of Chern 
number is straightforward and the quantized pumping 
approach in the topologically nontrivial (trivial) cases, 
corresponding to Sn x = 1 (Sn x = 0), is robust against 
weak perturbations under realistic conditions [561 [57]. 

Finally, we note that the probe of topological edge 
states in a cold atom system, such as those shown in 
Fig. 2 and Fig. 3, are also cumbersome. These edge 
states will generally be washed out by the smooth har¬ 
monic potential and thus one may not be able to distin¬ 
guish them from the bulk states. This problem could be 
circumvented potentially by using a steep confining po¬ 
tential and imaging the edge states from Bragg signals 
by means of specifically tuned Raman transitions m or 
from their dynamics after suddenly removing the poten¬ 
tial [62] . If the analogous Fermi line states in Fig. 3 can 
be created, one may further explore the induced topolog¬ 
ical current in the presence of an additional gauge field 
acting on the two Weyl points, which is used to break 
their central symmetry. The simplest example of such a 
gauge field can be generated by a constant Peierls phase 
for hopping along the z axis induced by the laser-atom 
coupling. 


V. CONCLUSIONS 

In summary, we have proposed an experimental scheme 
to simulate and explore 3D topological WSMs with cold 
atoms in a 2D square optical lattice subjected to an arti¬ 
ficial dimension from an external cyclical parameter. We 
have shown that this system is able to describe the essen¬ 
tial physics of WSMs with tunable Weyl points, and have 
investigated the relevant topological properties by calcu¬ 
lating the Chern number and the gapless edge states. 
Furthermore, we have proposed practical methods for 































the experimental detection of the mimicked Weyl points 
and the characteristic topological invariant in this cold 
atom system. Considering that all the ingredients to 
implement our scheme in the optical lattice have been 
achieved in the recent experiments, it is anticipated that 
the presented proposal will be tested in an experiment 
in the near future. Our proposed system would provide 
a promising platform for elaborating the intrinsic exotic 
physics of WSMs that are elusive in nature. 
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